Event-by-event simulation of nonclassical effects in two-photon interference experiments* 
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A corpuscular simulation model for second-order intensity interference phenomena is discussed. It is shown 
that both the visibility f = 1/2 predicted for two-photon interference experiments with two independent sources 
and the visibility 'f = 1 predicted for two-photon interference experiments with a parametric down-conversion 
source can be explained in terms of a locally causal, modular, adaptive, corpuscular, classical (non-Hamiltonian) 
dynamical system. Hence, there is no need to invoke quantum theory to explain the so-called nonclassical effects 
in the interference of signal and idler photons in parametric-down conversion. A revision of the commonly 
accepted criterion of the nonclassical nature of light is needed. 
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INTRODUCTION 

In classical optics, interference is known to be a phe- 
nomenon in which two waves are superimposed, resulting in 
a wave with bigger or smaller amplitude. Observed for the 
first time in Young's two-slit experiment in 1803 HI, it played 
an important role in the general acceptance of the wave char- 
acter of light. In quantum theory, interference in the two-slit 
experiment with electrons, large molecules, photons and other 
so-called quantum particles is considered to demonstrate the 
wave-particle duality of these quantum particles. In fact, ac- 
cording to Feynman the observation that the interference pat- 
tern in the two-slit experiment with electrons is built up de- 
tection event by detection event is a phenomenon which is 
"impossible, absolutely impossible, to explain in any classical 
way, and which has in it the heart of quantum mechanics" J2]. 
He referred to the interference of single electrons as "the only 
mystery" of quantum mechanics J2l. 

In general, a classical optical interference experiment con- 
sists of several classical light sources (not necessarily primary 
sources) and several detectors which measure the resulting 
light intensity at various positions. Adding equipment that 
accumulates the time average of the product of the detector 
signals allows for the measurement of the second and higher 
order intensity correlations. In quantum optics, the sources are 
replaced by single photon sources (the primary source com- 
monly said to create single photons or -photon entangled 
states with N >2) and single photon detectors. A coincidence 
circuit is added to the experimental setup to measure coinci- 
dences in the photon counts. 

In this paper we limit the number of sources and the number 
of detectors to two. Interference is then characterized by the 
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dependence of the resulting light intensity or of the second or- 
der intensity correlations on certain phase shifts. The Hanbury 
Brown-Twiss (HBT) effect was one of the first observations 
that demonstrated interference in the intensity-intensity corre- 
lation functions yfl. HBT showed that under conditions for 
which the usual two-beam interference fringes measured by 
each of the two detectors vanish, the correlated intensities of 
the two-detectors can still show interference fringes. For two 
completely independent sources, be it classical light sources 
or single photon sources, the visibility of this second-order in- 
tensity interference has an upperboundof 1/2 [4]. For primary 
sources producing correlated photon pairs, such as parametric 
down-converting sources, the two sources in an HBT-type of 
experiment can no longer be considered to be independent. In 
that case the two sources are considered to emit exactly one 
photon of the correlated pair simultaneously. Such sources 
provide a 100% visibility of the second-order intensity corre- 
lation, exceeding the 50% limit which is a commonly accepted 
criterion of nonclassicality [4]. The first experiment devoted 
to demonstrate nonclassical second-order intensity interfer- 
ence effects in the absence of first order intensity interfer- 
ence is probably the Ghosh-Mandel two-photon interference 
experiment of 1987 [5]. However, the effect is not limited to 
photons. Second-order intensity interference effects have also 
been observed in two-atom interference experiments |6|-|9|] in 
which an expanding cloud of cooled atoms acts as a source, 
multi-channel plate(s) detect the arrival and position of a par- 
ticle, and time-coincidence techniques are employed to obtain 
the two-particle correlations. Also in Hanbury Brown-Twiss 
type of experiments with electrons second-order intensity in- 
terference effects have been observed lioi 11]. As well in- 
tensity interference in the two-slit experiment as second-order 
intensity interference in Hanbury Brown and Twiss-type of 
experiments is attributed to the dual wave-particle character 
of the quantum particles. 
In previous work 

IBS 

we have demonstrated, using an 
event-based corpuscular model, that interference is not neces- 
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sarily a signature of the presence of waves of some kind but 
can also appear as the collective result of particles which at 
any time do not directly interact with each other. In general, 
the event-based approach deals with the fact that experiments 
yield definite results, such as for example the individual de- 
tector clicks that build up an interference pattern. We call 
these definite results "events". Instead of trying to fit the ex- 
istence of these events in some formal, mathematical theory, 
in the event-based approach the paradigm is changed by di- 
rectly searching for the rules that transform events into other 
events and, which by repeated application, yield frequency 
distributions of events that agree with those predicted by clas- 
sical wave or quantum theory. Obviously, such rules cannot 
be derived from quantum theory or, as a matter of fact, of any 
theory that is probabilistic in nature simply because these the- 
ories do not entail a procedure (= algorithm) to produce events 
themselves. 

The event-based approach has successfully been used to 
perform discrete-event simulations of the single beam splitter 
and Mach-Zehnder interferometer experiment of Grangier et 
al. 112 111 (see Refs. Il2l4l4ll) . Wheeler's delayed choice exper- 
iment of Jacques et aJ. 12211 (see Refs . il44l6ll ). the quantum 
eraser experiment of Schwindt et aJ. 12311 (see Ref. jl4ll7ll ). 
double-slit and two-beam single-photon interference experi- 
ments and the single-photon interference experiment with a 
Fresnel biprism of Jacques et aJ. [24] (see Ref. B14l 11810 . 
quantum cryptography protocols (see Ref. 12511 ). the Han- 
bury Brown-Twiss experiment of Agafonov et aJ. J26ll (see 
Ref. 11141 Il9ll ). universal quantum computation (see Ref. 1127 , 
2811 ). Einstein-Podolsky-Rosen-Bohm-type of experiments of 
Aspect et aJ. EH and Weihs et al (see Refs. JT1[32- 
3711 ). and the propagation of electromagnetic plane waves 
through homogeneous thin films and stratified media (see 
Ref. B14ll38ll ). An extensive review of the simulation method 
and its applications is given in Ref. [ 14]. Proposals for single- 
particle experiments to test specific aspects of the event-based 
approach are discussed in Refs. [ 1 



In this paper, we demonstrate that the second-order inten- 
sity interference with visibility 1/2 in a HBT experiment with 
two independent single photon sources and with visibility 1 
in the Ghosh-Mandel experiment can be entirely explained 
in terms of an event-based model, that is in terms of a lo- 
cally causal, modular, adaptive, classical (non-Hamiltonian) 
dynamical system. Hence, there is no need to invoke quantum 
theory to explain the observations and the commonly accepted 
criterion of the nonclassical nature of light needs to be revised. 

SECOND-ORDER INTENSITY INTERFERENCE 

In the context of the Ghosh-Mandel experiment, see Fig.Q] 
we may view the two mirrors as the two sources that produce 
two overlapping beams of photons. Hence, conceptually, this 
experiment can be simplified as shown in Fig. [2 which is the 
schematic diagram of a HBT experiment [4]. 

A HBT experiment is nothing but a two-beam experiment 
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FIG. 1: Diagram of the Ghosh-Mandel interference experiment |5]. 
A source emits pairs of single-photons through spontaneous down- 
conversion in a LiI03 crystal. These photons leave the source in 
different directions. Mirrors redirect the photons to the interference 
filter and a lens. The two beams overlap at a distance of about lm 
from the crystal. The resulting image is magnified by a lens and two 
movable glass pieces are used to collect and redirect the photons to 
the single-photon detectors Do and D\, the signals of which are fed 
into a coincidence counter CC. 
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FIG. 2: Schematic diagram of the Ghosh-Mandel experiment. Single 
photons emitted from point sources Sq and Si positioned at the y axis 
and separated by a center-to-center distance d are registered by two 
detectors Dq and D\ positioned on a line at a distance X from the 
y axis. The time of flight for each of the four possible paths from 
source S m to detector D„ is denoted by T,„,„ where m.n — 0, 1. 



with two independent sources and two detectors. The two 
sources are positioned along the y-axis and are separated by 
a center-to-center distance d. The two detectors are placed on 
a line at a distance X from the y-axis. Assume that source S m 
(m = 0, 1) emits coherent light of frequency / and produces a 
wave with amplitude A m e'^ m (A m and <j) m real). For simplicity 
of presentation, we assume that Aq = A\ = A. According to 
Maxwell's theory, the total wave amplitude B n on detector n 



B n =A(e 



(1) 



where the time of flight for each of the four possible paths 
from source S m to detector D n is denoted by T mM where m,n = 
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0,1. The light intensity /„ = \B„ 2 on detector D„ is given by 

/„ = 2A 2 { 1 + cos [fo - 0! + 27T/(r ,„ - r lin )] } . (2) 

If the phase difference 0o — 01 in Eq. (|2]i is fixed, the usual 
two-beam (first-order) interference fringes are observed. 

The essence of the HBT experiment is that if the phase dif- 
ference 0o — 01 is a random variable (uniformly distributed 
over the interval [0,2jt[) as a function of observation time, 
these first-order interference fringes vanish because 



(/„}=2A 2 , 



(3) 



where (.) denotes the average over the variables 0o and <j>i. 
However, the average of the product of the intensities Iq and 
I\ is given by 



(M>=4A 4 fl + icos2*/ArY 



(4) 



where Ar = (7qq — T\ q) — (Tq \ — T\ \). Accordingly, the 
intensity-intensity correlation Eq. (01 exhibits second-order 
interference fringes, a manifestation of the so-called HBT ef- 
fect. From Eqs. (01 and it follows that the visibility of the 
signal /, defined by 



max(7) — min(7) 
max (7) +min(7) 



(5) 



is given by Y = and Y = 1/2 for the first-order and second- 
order intensity interference, respectively. 

Treating the electromagnetic field as a collection of bosons 
changes Eq. © into J4J] 



bosons 



:4A 4 (l+cos27z;/Ar). 



(6) 



Clearly, for bosons, the visibility of the second-order intensity 
interference is "f = 1 . 

Considering the situation in which the two independent 
sources So and S\ are replaced by sources that emit simultane- 
ously exactly one photon of a correlated photon pair emitted 
by a parametric down-conversion source gives a similar ex- 
pression for the average of the product of the intensities Iq 
and I\ as given by Eq.(|6]l. Hence, also in this case Y = \ for 
the second-order intensity interference. 

In the two-beam experiment interference appears in its most 
pure form because the phenomenon of diffraction is absent. If 
we assume that the detectors cannot communicate with each 
other, that there is no direct communication between the par- 
ticles involved and that it is indeed true that individual pairs 
of particles build up the interference pattern one by one, just 
looking at Fig.|2]leads to the logically unescapable conclusion 
that the interference can only be due to the internal operation 
of the detector 13911 . Detectors that simply count the incoming 
photons are not sufficient to explain the appearance of an in- 
terference pattern and apart from the detectors there is nothing 
else that can cause the interference pattern to appear. We now 
discuss an event-based model of a detector that can cope with 
this problem IU4I1 . 



SIMULATION MODEL 

The model discussed in this paper builds on our earlier 
work Ill2l - ll4ll27ll . In short, in our simulation approach, a pho- 
ton is viewed as a messenger that carries a message and ma- 
terial is regarded as a message processor. Evidently, the mes- 
senger itself can be thought of as a particle. For the present 
purpose, it suffices to encode in the message, the time of flight 
of the particle. The interaction of the photons with material 
translates into a processing unit receiving, manipulating and 
sending out messages. Note that we explicitly prohibit two 
particles from communicating directly and that interference 
results from the processing of individual particles only 
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We now explicitly describe the model, that is we specify the 
message carried by the messengers, the algorithm for simulat- 
ing a detector ( = processing unit), and the simulation proce- 
dure itself. 

Messenger: The messenger can be regarded as a parti- 
cle which travels with velocity c in the direction q/q. Each 
messenger carries with it a harmonic oscillator which vibrates 
with frequency /. It may be tempting to view the messenger 
with its message as a plane wave with wave vector q, the os- 
cillator being one of the two electric field components in the 
plane orthogonal to q. However, this analogy is superfluous 
and should not be stretched too far. As there is no commu- 
nication/interaction between the messengers there is no wave 
equation (i.e. no partial differential equation) that enforces a 
relation between the messages carried by different messen- 
gers. Indeed, the oscillator carried by a messenger never inter- 
acts with the oscillator of another messenger, hence the mo- 
tion of these pairs of oscillators is not governed by a wave 
equation. Naively, one might imagine the oscillators tracing 
out a wavy pattern as they travel through space. However, 
as there is no relation between the times at which the mes- 
sengers leave the source, it is impossible to characterize all 
these traces by a field that depends on one set of space-time 
coordinates, as required for a wave theory. It is convenient 
(though not essential) to represent the message, that is the os- 
cillator, by a two-dimensional unit vector y = (cos y, sin y/) 
where \j/ = 2% ft + 8. Here, t is the time of flight of the parti- 
cle and 8 is a phase shift. Pictorially, the message is nothing 
but a representation of the hand of a clock which rotates with 
period 1 // and is running ahead by a time related to the phase 
8. A processing unit has access to this data and may use the 
messenger's internal clock to determine how long it took for 
the messenger to reach the unit. 

Source: A source creates a messenger with its phase 8 set 
to some randomly chosen value. Initially its time of flight t 
is zero as it is determined by the arrival of the messenger at 
a processing unit. A pseudo-random number determines to 
which detector the messenger travels. 

Single-photon detector: In reality, photon detection is the 
result of a complicated interplay of different physical pro- 
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FIG. 3: Diagram of the event-based detector model defined by 
Eqs. - d 10b - The detection unit consists of an input stage, which is 
a deterministic learning machine (DLM), a transformation stage and 
an output stage. The input stage has K input channels at which a mes- 
sage y, being a two-component vector, can arrive, K corresponding 
internal registers in which the incoming message can be stored 
and one internal A"-component vector x, responsible for the learning. 
The transformation stage generates a message T, a two-component 
vector, based on all information available in the input stage. The out- 
put stage takes the message T as input and generates an output signal 
z representing a "click" or "no click" on output channel or 1, re- 
spectively. The detection unit processes one message at a time. The 
solid lines indicate the input and output channels of the processing 
unit and the dashed lines indicate the data flow within the processing 
unit. 



In essence, a light detector consists of material that absorbs 
light. The electric charges that result from the absorption pro- 
cess are then amplified, chemically in the case of a photo- 
graphic plate or electronically in the case of photodiodes or 
photomultipliers. In the case of photomultipliers or photodi- 
odes, once a photon has been absorbed (and its energy "dis- 
sipated" in the detector material) an amplification mechanism 
(which requires external power/energy) generates an electric 
current (provided by an external current source) Il40l 14111 . The 
resulting signal is compared with a threshold that is set by 
the experimenter and the photon is said to have been detected 
if the signal exceeds this threshold I140L I41I1 . In the case of 
photographic plates, the chemical process that occurs when 
photons are absorbed and the subsequent chemical reactions 
that renders visible the image serve similar purposes. 

Photon detectors, such as a photographic plate of CCD ar- 
rays, consist of many identical detection units each having a 
predefined spatial window in which they can detect photons. 
In what follows, each of these identical detection units will 
be referred to as a detector. By construction, these detector 
units operate completely independently from and also do not 
communicate with each other. 

An event-based model for the detector cannot be "derived" 
from quantum theory simply because quantum theory has 
nothing to say about individual events but predicts the fre- 
quencies of their observation only B42I1 . Therefore, any model 
for the detector that operates on the level of single events must 
necessarily appear as "ad hoc" from the viewpoint of quantum 
theory. The event-based detector model that we employ in this 



paper should not be regarded as a realistic model for say, a 
photomultiplier or a photographic plate and the chemical pro- 
cess that renders the image. In the spirit of Occam's razor, the 
very simple event-based model captures the salient features of 
ideal (i.e. 100% efficient) single-photon detectors. 

The key element of the event-by-event approach is a pro- 
cessing unit that is adaptive, that is it can learn from the 
messengers that arrive at its input ports Il2l4l4il . The dia- 
gram of an event-based detection unit is depicted in Fig. [3] 
It consists of an input stage called deterministic learning ma- 
chine (DLM) BI2II13II . a transformation stage, and an output 
stage. The processing unit should act as a detector for in- 
dividual messengers which may come from several different 
directions. Therefore, as can be seen from the schematic dia- 
gram depicted in Fig. [3]this processing unit has K input ports, 
a parameter that allows the machine to resolve K different di- 
rections. 

Input stage: Representing the arrival of a messenger at port 
1 < k < K by the vector v = (vi, . . . , vk) t with v; = 5,- \ (i = 
l,...K) the internal vector is updated according to the rule 



yx+(l-y)v, 



(7) 



where x = (x\, . . . ,Xk) T , Lf=i x k = 1> an d < 7 < 1. The 
elements of the incoming message y are written in internal 
register 



(8) 



while all the other Y, (2 ^ k) registers remain unchanged. 
Thus, each time a messenger arrives at one of the input ports, 
say k, the DLM updates all the elements of the internal vector 
x, overwrites the data in the register Yj- while the content of 
all other Y registers remains the same. 

Transformation stage: The output message generated by 
the transformation stage is 



T = x Y 



K 



(9) 



which is a two-component vector. Note that |T| < 1. 

Output stage: As in all previous event-based models for the 
optical components, the output stage generates a binary output 
signal z = 0, 1 but the output message does not represent a 
photon: It represents a "no click" or "click" if z — or z = 1, 
respectively. To implement this functionality, we define 



0(|Tf 



(10) 



where ©(.) is the unit step function and < < 1 are uniform 
pseudo-random numbers (which are different for each event). 
The parameter < 7 < 1 can be used to control the operational 
mode of the unit. From Eq. ( fTOb it follows that the frequency 
of z — 1 events depends on the length of the internal vector T. 

Note that in contrast to experiment, in a simulation, we 
could register both the z — and z = 1 events. Then the sum 
of the z = and z = 1 events is equal to the number of input 
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messages. In real experiments, only z = 1 events are taken 
as evidence that a photon has been detected. Therefore, we 
define the total detector count by 



N 

Account = Zl , 
1=1 



(ID 



where N is the number of messages received and I labels the 
events. In other words, Amount is the total number of one's 
generated by the detector unit. 

Comparing the number of ad hoc assumptions and un- 
known functions that enter quantum theoretical treatments of 
photon detectors [41] with the two parameters y and K of the 
event-based detector model, the latter has the virtue of being 
extremely simple while providing a description of the detec- 
tion process at the level of detail, the single events, which in 
any case is outside the scope of quantum theory. 

Simulation procedure: Before the simulation starts we set 
x = (1,0, . . . ,0) r and we use pseudo-random numbers Si to 
set Y* = (cos2^^',sin2^^') for k = 1, . . . ,K. Next, we gen- 
erate iVtot pairs of messengers, send them to the detectors, de- 
termine the detector count Amount a t Do and D\ and count the 
coincidences. In the simulation always two messengers travel 
to the detectors, one generated at source Sq and one at source 
Si . Hence, once a pair of messengers is generated a detector 
can generate no click, one click or two clicks. Only when both 
detectors generate a click the coincidence count Moincindence is 
enhanced by one. 



SIMULATION RESULTS 
Detection efficiency 

The efficiency of the detector model is determined by sim- 
ulating an experiment that measures the detector efficiency, 
which for a single-photon detector is defined as the overall 
probability of registering a count if a photon arrives at the de- 
tector [40]. In such an experiment a point source emitting 
single particles is placed far away from a single detector. As 
all particles that reach the detector have the same time of flight 
(to a very good approximation), all the particles that arrive at 
the detector will carry nearly the same message y which is en- 
coding the time of flight. Furthermore, they arrive at the same 
input port, say q. As a result x (see Eq. (0) rapidly converges 
to the vector with X{ —> 8i. q and, as y is a unit vector, we have 
|T| 1, implying that the detector clicks almost every time 
a photon arrives. Thus, for our detector model, the detection 
efficiency as defined for real detectors [40] is very close to 
100% (results not shown). 



Hanbury Brown-Twiss experiment 

In Fig. |4]we present the simulation results for the HBT ex- 
periment depicted in Fig. [2] For simplicity, we have put detec- 
tor Dq at (X,Q) and plot the single detector and coincidence 
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FIG. 4: Simulation data of the single-particle and two-particle counts 
for the HBT experiment depicted in Fig. [2] Red open circles (Blue 
open triangles): results for the counts Amount of detector Dq (Di), 
showing that there is no second-order intensity interference. Red 
closed circles: results for the coincidence counts A^ i nc i(j ence . The 
dashed and solid lines represent the theoretical predictions N tot /2 
and Eq. J 1 3b for the single detector and coincidence counts, respec- 
tively. Simulation parameters : A^t = 2 x 1 6 events per y\f j c- value, 
N F = 50, X = 100000c//, d = 2000c//, y = 0.99 and K = 2. 



counts as a function of the y-position of detector D\, In each 
simulation step, both sources 5b and Si create a messenger 
with some randomly chosen phase being the only initial con- 
tent of the messages y m (m = 0, 1). The phases are kept fixed 
for Np successive pairs of messengers. The total number of 
emitted pairs is denoted by N tot . Two pseudo-random num- 
bers are used to determine whether the messengers travel to 
detector Do or D\. The time of flight for the messenger trav- 
elling from source S m to detector D„ is given by 



T — 

1 m,n 



^/X 2 + ((l-2m)d/2-y n y 



(12) 



where m,n — 0, 1. The time of flight T Mt „ is added to the mes- 
sage y m before the message is processed by the corresponding 
detector D„. The messages are the only input to the event- 
based model. As Fig.|4]shows, averaging over the randomness 
in the initial messages (random phases) wipes out all interfer- 
ence fringes in the single-detector counts, in agreement with 
Maxwell's theory. We find that the number of single-detector 
counts Amount fluctuates around Mot/2, as expected from wave 
theory. Similarly, the data for the coincidence counts are in ex- 
cellent agreement with the theoretical prediction for the sim- 
ulation model 

Mot /. 1 



N, 



coincidence 



: cos2nfAT 



(13) 



and, disregarding the prefactor Mot/8, also in qualitative 
agreement with the predictions of wave theory. 

For simplicity, we have confined the above presentation to 
the case of a definite polarization. Simulations with randomly 
varying polarization (results not shown) are also in concert 
with Maxwell's theory. 
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Ghosh-Mandel experiment 

From Eq. ( fl"3l >. it follows that the visibility of the interfer- 
ence fringes, defined by 



max (N coincidence ) min(iV coincidence ) 
max (A^ co i nc idence ) H~ ^i^(^coincidence) 



(14) 



cannot exceed 1/2. It seems commonly accepted that the vis- 
ibility of a two-photon interference experiment exceeding 1/2 
is a signature of the nonclassical nature of light. 

As two-photon interference experiments, such as the Gosh- 
Mandel experiment J^t], employ time-coincidence to measure 
the intensity-intensity correlations, it is quite natural to expect 
that a model that purports to explain the observations accounts 
for the time delay that occurs between the time at which a par- 
ticle arrives at a detector and the actual click of that detector. 
In quantum theory, time is not an observable and can there- 
fore not be computed within the theory proper. Hence there 
is no way that these time delays, which are being measured, 
can be accounted for by quantum theory. Consequently, any 
phenomenon that depends on these time delays must find an 
explanation outside the realm of quantum theory (as it is for- 
mulated to date). 

It is straightforward to add a time-delay mechanism to the 
event-based model of the detector. For simplicity, let us as- 
sume that the time delay for the detector click is given by 



^delay 



- Wl-|T| 2 ) /J ln^ 



(15) 



where < M < 1 is a pseudo-random number, and T is given 
by Eq. ([9j. The time scale T max and the exponent h are free 
parameters of the time-delay model. Note that fdeiay — %n,n is a 
pseudo-random variable drawn from an exponential distribu- 
tion with mean r max (l — |T| 2 ) ft . Coincidences are counted by 
comparing the difference between the delay times of detectors 
Do and D\ with a time window W. 

From the simulation results presented in Fig. [5J it is clear 
that by taking into account that there are fluctuations in the 
time delay that depend on the time of flight and the inter- 
nal state of the detector, the visibility changes from — 1/2 
to "V » 1 . The simulation data is represented (very) well by 
Mot/2 and 



N' 

1 "count 



^coincidence ' 



: a\N t< 



(1 + cos27r/Ar), 



(16) 



where the prime indicates that the model incoporates the time- 
delay mechanism and a' 4 is a fitting parameter which depends 
on the details of the time-delay mechanism. As expected, the 
use of a narrow time window leads to a signicant reduction 
(by a factor a' 4 = 0.077) of the total coincidence count. These 
results demonstrates that a purely classical corpuscular model 
of a two-photon interference experiment can yield visibilities 
that are close to one. Hence, the commonly accepted criterion 
of the nonclassical nature of light needs to be revised. 

The time delay model Eq. ( TT3T > is perhaps one of the sim- 
plest that yield interesting results but it is by no means unique 
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FIG. 5: Simulation data of the single-particle and two-particle counts 
for the HBT experiment depicted in Fig. [2] generated by the same 
event-based that produced the data of Fig.[4]extended with the time- 
delay model Eq. dl5t . Simulation parameters: r max // = 1000, 
W/f = 1, h = 8. The dashed and solid lines are least-square fits 
to a'2N tot and a' 4 N tot (l + b' 4 cos2KfAt) for the single detector and 
coincidence counts, Amount and A' co ; nc i ( ience> respectively. The values 
of the fitting parameters are a' 2 = 0.502, a\ = 0.077 and b' 4 = 0.974. 



and can only be scrutinized on the basis of accurate experi- 
mental data which, unfortunately, do not seem to be available 
thus far. 



Bosons 

If we exclude the possibility that the two sources send their 
particles to the same detector, the event-based approach pro- 
duces results that are reminiscent of the quantum theoretical 
description in terms of bosons. In Fig. [6] we present the re- 
sults of such a simulation, using the same model parameters 
as those used to produce the results of Fig. [5] From Figs.|5]and 
[6] it is clear that the maximum amplitude of the two-particle 
interference signal of the latter is two times larger than that of 
the former (the "classical" case), as expected for bosons. The 



simulation data is represented (very) well by 



n: 



if 

coincidence 



(1 + cos2?r/Ar), 



-Mot/2 
(17) 



where the double prime indicates that the model incoporates 
the time-delay mechanism and that the possibility that the two 
sources send their particles to the same detector has been ex- 
cluded. 



Non-monochromatic sources 

All the results presented above have been obtained by as- 
suming that the beams of particles are strictly monochromatic, 
meaning that the frequency / of the oscillators carried by the 
particles is fixed. A more realistic simulation of the pairs of 
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the source. As well the interference experiment with a classi- 
cal light source as the interference experiment with the para- 
metric down-conversion source can be explained entirely in 
terms of a classical corpuscular model. 

Elsewhere, we have shown that third order intensity inter- 
ference in a Hanbury Brown-Twiss type of experiment with 
two sources emitting uncorrected single photons can be mod- 
eled by an event-based model as well I19B . Simulation of an 
interference experiment with a three-photon source is left for 
future research. 



FIG. 6: Same as Fig. [5] except that the two sources never send their 
particles to the same detector, mimicking bosons (see text). The 
values of the fitting parameters are a" = 0.502, a'^ = 0.154 and 



b'l = 0.985. 



photons created by the parametric down-conversion process 
requires that the frequencies fa and fa of the messages car- 
ried by the pair of particles satisfy energy conservation, mean- 
ing that h + fi = fo where /o is the frequency of the pump 



beam 11411 143M45I1 . It is straightforward to draw the frequen- 
cies fi (and therefore fa = fn — fa ) from a specified distribu- 



tion, such as a Lorentzian H41 , 



In the simulation, each 



created particle pair would then correspond to one message 
characterized by a frequency fa and another one by frequency 
fa. The detectors simply sum all the contributions (taking into 
account the differences in the factors f m T m ^ n ), resulting in a 
reduction of the visibility, just as in the wave mechanical pic- 
ture. 



CONCLUSION 

We have shown that the so-called nonclassical effects ob- 
served in two-photon interference experiments with a para- 
metric down-conversion source can be explained in terms of a 
locally causal, modular, adaptive, corpuscular, classical (non- 
Hamiltonian) dynamical system. The high visibility, "V = 1, 
in this type of experiment is commonly considered as a signa- 
ture of two-photon light, in contrast to the visibility "f = 1/2 
obtained in a similar experiment with a classical light source. 
On the other hand, according to Ref. B26I1 . the existence of 
high-visibility interference in the third and higher orders in 
the intensity cannot be considered as a signature of three- or 
four-photon interference, because high-visibility interference 
is also observed in Hanbury Brown-Twiss type interference 
experiments with classical light. Hence, although the case 
of second-order intensity interference seemed to be different 
from the higher orders, we have demonstrated that also for the 
second order intensity interference the value of the visibility 
cannot be used to say anything about the quantum character of 
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